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Abstract 



5 : We show that the low frequency ' long wavelen s th d ^ namics of the P hase of the P air field 

for a BCS-type s-wave superconductor at T=0 is equivalent to that of a time-dependent 



non- linear Schrodinger Lagrangian (TDNLSL), when terms required by Galilean invariance 
are included. If the modulus of the pair field is also allowed to vary, the system is equivalent 
to two coupled TDNLSL's. 
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The classic Ginzburg-Landau (GL) theory [1] is very successful [2, 3] in describing a large 
class of static superconducting phenomena near the critical temperature T c , and its form 
was established by Gorkov [4] shortly after the microscopic BCS theory [5]. Subsequently, a 
number of attempts [2,3,6,7] were made to obtain a generalized GL theory for time-dependent 
phenomena, and for temperatures well below T c , but a consensus has still not been reached 
on the form of such a theory at T=0. In this paper we shall show that the low frequencylong 
wavelength dynamics of the phase of the pair field for a BCS-type s-wave superconductor 
at T=0 is equivalent to a time-dependent non-linear Schrodinger Lagrangian (TDNLSL). 
At first sight, this result might seem almost obvious: after all, the energy density in GL 
theory looks formally like that of a non-linear Schrodinger theory so that it seems natural to 
extend it to the corresponding time-dependent theory as, indeed, Feynman assumed [8] in his 
discussion of the dynamics of superconductors and of the Josephson effects. Yet neither the 
earlier discussions [2, 3, 6], nor recent work based on the effective action formalism of quantum 
field theory [7,9], appears to lead to this conclusion. This is in contrast to the case of a Bose 
superfluid, such as 4 He, which is well described by a TDNLSL near T=0[10]. Indeed, there is 
considerable current interest in probing the relationship and "crossover" between BCS and 
Bose superfluidity [11]. Our result implies that both are fundamentally the same, at least 
near T=0 in the clean limit ; in particular, the existence of the Magnus force for a vortex 
line in a superconductor follows naturally. The last point is pertinent to the discussion of 
vortex dynamics in superconductors within the effective theory formulation [12]. 

Three of the present authors have, in fact, recently shown [13] that the motion of the 
condensate is described by a non-linear Schrodinger equation at T=0, using a density matrix 
approach and the Born-Oppenheimer approximation. But this left open the question how 
this could be reconciled with the earlier work [2,3,6,7,9], which was generally based on field 
theory (or Green function) techniques, and apparently led to a quite different result. The 



solution of this problem is contained in the present paper, and it is essentially very simple. 
In [13] the further approximation was made that the modulus of the energy gap (or pair 
field) is constant. If this approximation is made in the field theory approach, one can 
derive [9] from BCS theory an effective action for the phase 9(x) of the pair field (i.e. the 
Goldstone mode), which one then expands up to a certain order of derivatives. The resulting 
Lagrangian L e ff(9) is the same as that previously proposed on symmetry grounds [14], and 
also corresponds precisely to the early results of Kemoklidze and Pitaevskii [15], who started 
from Gorkov's equations [5]. We shall show that the dynamics of 9(x) as given by L e ff(9) 
can be re- written in terms of a TDNLSL, which is equivalent to a particular case of the 
general non-linear Schrodinger theory derived in [13] under the same approximation. 

We also extend this to include variations in the modulus of the pair fields, and show 
that the dynamics is then that of two coupled TDNLSL 's. The thread that unites all these 
approaches is ultimately Galilean invariance. Since the microscopic starting point is always 
Galilean invariant, one expects any theory based on an effective action to preserve this 
symmetry, a point emphasized in Ref. 15, and the Schrodinger Lagrangian is the simplest 
such available. 

We begin with the BCS Lagrangian for s-wave pairing and in the absence of external 
fields: 

L = J2^*A X ) Udt + — + /2J ip a (x) +# T *(x)^J(x)^|(x)^ T (a;) (1) 

where ip a describes electrons with spin a = (|, j), ^ — k 2 F /2m is the Fermi energy in the 
normal state, and x = (x, t). Introducing the auxiliary (pair) fields A(x) and A*(x), and 
integrating out the electron fields, one obtains the effective action 

S[A, A*] = -%Tr In G' 1 - - f d 4 x\ A(x) | 2 (2) 



where the Nambu Green function satisfies 

( d A(x) ^ 



G(x,x') = 5(x-x') (3) 



\ A*(x) 2 j 

with 0\ = id t + V 2 /2m + //, 2 — id t — V 2 /2m — fi, and Tr includes interval and space- 
time indices. To obtain from (2) an effective Lagrangian in terms of the degrees of freedom 
represented by A, a possible procedure [7] is to set A(x) = Ao + A'(x) where Ao is the 
position of the minimum of S for space-time independent A, and where A' is assumed to be 
slowly varying in both space and time. One then expands Tr In G" 1 in powers of derivatives 
of A'. There are, however, two (related) objections to this. First, we are dealing with the 
spontaneous breaking of a local U(l) phase invariance, which implies that at a temperature 
far from the transition temperature, the most important degree of freedom is the phase of A, 
which is the relevant Goldstone field. It is this field, rather than the real and/or imaginary 
parts of A, which should carry the low frequency and long wavelength dynamics. Secondly, 
the ansatz A(x) = A + A'(x) violates the Galilean invariance possessed by (1), which implies 
[15] that 

A(r — vt, t) exp(2imv • r — imv 2 t) (4) 

should satisfy the same equation of motion as A(r, t). We shall return to the question of 
Galilean invariance below. 
We therefore set 

A(x) = e w{x) \A(x)\ (5) 

and |A(x)| = | A | +5 \ A (x)\, where we are interested in the low frequency and long wavelength 
fluctuations of 6(x) and 5|A(x)|. However, although 5|A(x)|/|Ao| is expected to be small, 
and a simple expansion of the sort mentioned above for Tr In G~ l could easily be set up 
in terms of derivatives of 5|A(x)| if 9(x) were zero, it is crucial to recognize that 6(x) is 



not small in general, so that the phase factor in (5) cannot be expanded, but must be 
treated as a whole. This prevents a straightforward expansion of Tr In G~ x when (5) is 
substituted into (2). Fortunately, this difficulty can be easily circumvented [9,16]. Defining 
U(x) = exp(i9(x)r 3 /2), we can write 

Tr In G- 1 = Tr In G^UU' 1 = Tr In UG^U' 1 = Tr In G' 1 (6) 

where 



G- 1 = Gn Hi - G n E) 

(7) 
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and 

E = ~iV 2 9/Am - iV9 ■ V/2m + (6/2 + (V9) 2 /8m)r 3 - 5\A\ n. (8) 

Minimizing (2) with 9 = S\A\ = yields the usual gap equation for |A |. The dynamics of 
9 and IAI is contained in 
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S eff [9, S\A\] = iTr £ -(G S)" - - y | A\ 2 d 4 x , (9) 



where we note that S contains just 5|A(x)| and derivatives of 9(x), in terms of which (as- 
sumed small) quantities a useful expansion can be conducted, following standard techniques 
[17]. 

We now concentrate on 9(x), and set S\A\ = for the time being. Carrying out the 
calculation to the indicated order in derivatives we obtain 

L eff(0) = ~Po(0/2 + (W) 2 /8m) + JV(O)(0/2 + (W) 2 /8m) 2 (10) 

where po — kp/3n 2 is the electron density at T=0, N(0) is the density of states ( for one 
spin projection) at the Fermi surface and we have adopted a convenient normalization; note 



that iV(O) = po/2mVa where v a = vf/V3 is the velocity of propagation of the Bogoliubov- 
Anderson mode. Eqn. (10) is the same as Eqn. (4) in Ref.9 if 9 is replaced by 20 (see also 
Eqn.(2.6) of ref.14). The equation of motion for 9 which follows from (10) is 

| + v.j = o (ID 

where 

p = Po-N(0)(6 + (V6) 2 /4m) (12) 

and 

j = pV6/2m. (13) 

Equations (11)-(13) are, in fact, precisely those obtained (to this order in derivatives) by 
putting S\A\ = in Eqns.(21)-(23) of Ref.15. We now show how the dynamics contained in 
(11)-(13) can be reinterpreted in terms of a TDNLSL. 

We first remark that the forms of (11) and (13) strongly suggest that the quantities p and 
j have the physical significance of a number density and of a number current density, respec- 
tively In fact, simple linear response theory (assuming, as always, a derivative expansion) 
gives [18] 

5p « -N(0)9, j « p V9/2m (14) 

where Sp is the departure of the density from the equilibrium value p , and j is the number 
density current. We now consider the implications of Galilean invariance. For densities p 
and j obeying (11), Galilean symmetry requires 

p'(r', t') = p(r, t), j'(r', t') = j(r, t) - vp(r, t) (15) 

where r' = r—\t,t' = t; that is, p and j transform as t and r respectively. From the discussion 
around Eqn. (4) above, the behaviour of 9 under Galilean transformations is given by 

9'(r', t') = 0(r, t) + mv 2 t - 2mv • r, (16) 



from which it follows that 5p as given by (14) is not invariant, and that j transforms in- 
correctly. However, although 9 is not a Galilean invariant, the combination 9 + (V9) 2 /4m 
is. (This is actually the justification for singling out the terms given in (10) from the total 
n = 2 contribution in (9).) Thus the requirement that Sp be Galilean invariant leads pre- 
cisely to the expression (12) for p — p , which can now be identified as 5p. Similarly, if 
we replace po m the expression (14) for j by p, we find that j transforms correctly and the 
expression for j is that in (13). 

We are therefore led to seek a theory involving two fields p and 9, such that the equation 
of motion for 9 is (11) and that for p is (12), with p, 9 and j related by (13). Consider the 
TDNLSL 

L , = j^*?fj - -^V^* • W - V (17) 

4m 

where the mass has been chosen to be 2m, and V will be assumed to be a function of \ijj\ 
only. Let us set 

if> = v/p exp(^). (18) 

Then inserting (18) into (17) and discarding a total derivative, L^ becomes 

L f = - P 9 - p(V9) 2 /4m - (Vp) 2 /16mp - V{p). (19) 

The equation of motion for 9 is then (11), with j given by (13), while that for p is 

^ = _(0 + (W) 2 /4m) - (Vp) 2 /16mp 2 + V 2 p/8mp 2 . (20) 

dp 

If we now choose 

V = (p- p ) 2 /2N(0) (21) 

and solve (20) by expanding in derivatives, we recover precisely (12) at the relevant order. 

Thus the introduction of the auxiliary variable p, which can be expressed in terms of 9 
via its equation of motion, has allowed us to rewrite the dynamics of the Goldstone field 9, 
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as given by L e ff(0), in terms of the TDNLSL (17). The variable p is interpreted physically 
as the number density It must be stressed, however, that the wavefunction if) introduced in 
this way (see (18)) is quite distinct from the pair field A, despite the fact that they have the 
same phase 9. In our development thus far, the modulus |A| has been held fixed, whereas p 
varies; there is no simple relation between |A| and ^fp. 

The dynamical theory described by (17) is a special case of the general time-dependent 
non-linear Schrodinger theory for the condensate wavefunction derived in [13], in which the 
form of the potential V was not explicitly calculated. Here, it has been necessary to fix V 
in order to carry out the elimination of the variable p, to the required order in derivatives. 
(We remark that in [13] if) is normalised to the density of Cooper pairs, rather than, as here, 
to the electron density.) 

Before discussing the modifications caused by the inclusion of the field <5|A(x)|, we make 
one further comment on the Lagrangian (10). A simple alternative route to (10) is to start 
from a Lagrangian which describes just the Bogoliubov- Anderson mode, viz. 

L a = \& - \v 2 a {V9f. (22) 

Now (22) is clearly not invariant under (16). But, as we have seen, the combination 
6+ (V6») 2 /4m is invariant. Hence if in (22) we replace 9 2 by (9+ (V6») 2 /4m) 2 and (V6>) 2 by 
4m9 + (V9) 2 we will have a Galilean invariant Lagrangian; and the result of these replace- 
ments is just proportional to L e ff(9). Actually, the (V#) 2 term in (22) is of course invariant 
by itself, up to constants and a total derivative. Indeed, the term 9 introduced above, and 
present in L e ff(9), is also a total derivative and does not affect the equations of motion. 
Nevertheless it is important physically, as it ensures that the density p has the equilibrium 
value po [14]. 

In view of its relative unfamiliarity, it maybe worth noting that L e ff(9) (or equivalently 



Lty) embodies the usual phenomenology of superfluid dynamics at T=0 (see, for example 
[19,20]). We identify V#/2m with the superfluid velocity v s , and multiply p and j of (12), 
(13) by m to convert them to mass density and flux, p m and j m . Eq. (11) is then the law of 
mass conservation, following from the fact that L e ff does not depend explicitly on 6. Eq.(12) 
is equivalent [15] to Bernoulli's equation, if we make use of 5p ~ 2N(0)Sp and Sp ~ PoSp. 
Since L e ff does not depend explicitly on t, we have the energy conservation relation 

dE 

— + V • Q = (23) 

where using the canonical definitions (with suitable normalization), one finds 

E « -p m v 2 s , Q = j m (-v 2 + <5/i) , (24) 

and we have dropped a quantity of order Sp 5p in E. Finally, since L e ff is translation 
invariant we have the momentum conservation relation 



% + V • n = (25) 

at 

where the momentum flux density tensor is 

% = PmV s iV sj + Sp Sij. (26) 

Eq.(25) is equivalent to Euler's equation. In Ref.14, the proportionality between the mo- 
mentum density and the number current j (defined by dL/d(V6)), which is included in (25), 
was taken as a constraint on possible Lagrangians L. If L is a function solely of the Galilean 
scalar g = 9 + (V^) 2 /4m, then 

dL _dL dL dL V6 

~d9 ~~dg' d(W) ~ ~dg2m' ( ' 

Since 6 is a phase variable, we can interpret dL/d6 and dL/d(V6) as being proportional 
to a conserved number density p and number current density j respectively, so that (27) 
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becomes just (13). The momentum density is then automatically proportional to j. Once 
again, Galilean invariance is the essential principle. 

We now turn to the inclusion of the field 5|A(rr)|. L e ff(6, S\A\) can be extracted from 
(9), up to a given order in derivatives, but calculations rapidly become laborious. For our 
present purpose, we will simply use the result of Ref.15 which, using the normalization of 
(10), gives 



J eff\ 



-Po 



9/2 + 



(W) 2 , (Ve) 



+ N(0) 



im 



e/2 + 



im 



+ iV(0) 



9/2 + 



(W) 2 , (Ve) 



5ffl 



im 



Ve-Vfl 

4m 



2 2 



3iV(0)|A | 2 e 



(28) 



where e(x) = 5|A(x)]/(v3|Ao|) and we have retained corresponding terms in e and 9. The 

quadratic terms in e yield the amplitude collective mode with a finite gap found in Ref.18 

(and are also in agreement with the result of Ref.7); we have omitted higher powers of e. 

The term in e/2 is made Galilean invariant by the addition of Ve • V#/4m, since e(x) is a 

scalar. We now find that in order to rewrite (28) in "Schrodinger" form we have to introduce 

a second auxiliary density feld, which we call p e , in addition to the earlier density p, which 

now becomes pg. Thus four fields (9, e, pg and p e ) are required, and (28) is actually equivalent 

to two coupled Schrodinger equations. That is, the equations of motion for 9 and e which 

follow from (28) are identical to those arising from 

1 



£l/>l,</>2 = #1^1 



4m 



V^-V^i-(|Vir-po/2)/iV(0) 



+ # 2 >2 - t^W> 2 * • W> 2 - (|^ 2 | 2 - Po/2)/N(0) + ^iV(0)|A | 2 [/mln^iM)] 2 
4m " 2 



where 



V>i = ^{ Pe + p e )/2 exp(i(9 + e)), ^ 2 = y/(p e - p e )/2 exp(i(9 - e)). 
For example, corresponding to (12), we have 



„ = »- m ( i+ 5?!r + w 



4m 



4m 



09 



(29) 



(30) 



(31) 
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and 



*--W« + ^-^ 



Eqn.(29) represents a system of two TDNLSL's coupled via the "mass" term in (28). Ex- 
pressions for all the conserved quantities can be found as before, and will include quantum 
corrections to the semiclassical results of (23)- (26). 

The inclusion of electromagnetism in the above formalism is straightforward. Consider 
the formulation in terms of L e ff(6,e). Since 9 is the phase of a field with charge — 2e, 
gauge invariance implies that 9 and V# must appear in the combinations 9 — 2eA and 
V# + 2eA (e > 0), where A and A are the electromagnetic potentials. The field e, on 
the other hand, is electromagnetically neutral. The leading order electromagnetic charge 
and current densities are obtained by multiplying Sp and j in (22) by — e and making the 
above replacements for 9 and V#. One then obtains the usual results [7]. In terms of the 
Schrodinger formulation, one simply makes the expected minimal coupling substitutions: 
id t — » id t + 2eA and — iV — > — iV + 2eA in (14) or (29) (note from (30) that both ^ and 
■02 have charge — 2e). 

When the above analysis is extended to higher order derivative terms, it is clear on 
dimensional grounds that some characteristic scale must enter. In fact, such higher terms 
enter in the form d t /\A \ and i^V/|A | ~ £V (see for example Eqn.(35) of [9]), where £ is 
the coherence length. The basis of the expansion is therefore the usual assumption [15] that 
the characteristic frequency ou, and wavenumber k, of variations of A(x) satisfy uo <^i |Ao|, 
k <^ £ _1 . Indeed, (28) already yields a static solution for e which decays exponentially over 
a characteristic distance £/6. Such a solution is of the type expected far from a vortex core. 
Inclusion of appropriate higher derivative terms should make possible some predictions about 
the vortex core structure. 

The TDNLSL formulation provides, we believe, a simple and unifying framework for the 
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discussion of dynamical effects in BCS superconductors at T=0. Results which have been 
known for many years [2,3,6,15], as well as those obtained by quite different methods only 
recently [9,14], are all seen to be in agreement with each other, and with the TDNLSL 
formulation. 

Note Added: After completion of this work we received a copy of a preprint by Michael 
Stone [21], in which a similar conclusion is reached concerning the effective TDNLSL when 
|A| is constant. 
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